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INTRODUCTION 


BACKGROUND 

The  Military  Specification  for  Flying  Qualities  of  Piloted  Airplanes,  MIL-F-8785B,  (reference 
(a)),  was  developed  largely  from  flight  tests  of  classically  responding  unaugmented  aircraft.  Its 
quantitative  requirements  are  generally  expressed  in  terms  of  modal  approximations  which  can  be 
described  mathematically  by  first  or  second  order  linear  expressions.  Advancements  in  aerody¬ 
namics  and  complicated  control  system  augmentation  schemes,  prevalent  in  modern  aircraft  designs, 
have  resulted  in  responses  which  are  described  by  high  order  functions. 

In  an  attempt  to  utilize  the  existing  requirements  in  analyzing  advanced  aircraft/control  sys¬ 
tem  configurations,  the  concept  of  equivalent  systems,  has  been  introduced  (reference  (b) ) .  A 
digital  frequency  domain  equivalent  system  matching  technique  has  been  developed  by  Hodgkinson, 
et.  al.,  and  applied  to  the  high  order  representations  of  experimental  aircraft  (reference  (c)).  The 
approach  used  was  to  approximate  the  high  order  response  to  pilot  control  input  transfer  functions 
of  the  subject  aircraft  with  classical  low  order  transfer  functions  describing  the  specification  require¬ 
ments,  augmented  with  a  time  delay.  This  equivalent  time  delay  approximates  the  phase  lag  intro¬ 
duced  by  the  high  frequency  control  system  components.  Within  the  scope  of  the  initial  investiga¬ 
tions,  it  was  determined  that  the  linear  modal  requirements  of  MIL-F-8785B,  when  augmented  by 
a  requirement  on  time  delay,  are  appropriate  for  specifying  the  handling  qualities  of  the  advanced 
high  order  configurations  of  tomorrow's  airplanes  (reference  (d) ) .  This  approach  has  been  incor¬ 
porated  in  the  latest  revision  to  the  MIL-SPEC,  MIL-F-8785C,  and  the  proposed  MIL  standard  and 
Handbook,  reference  (e)  and  (f)  respectively,  which  states: 

"The  contractor  shall  define  equivalent  classical  systems  which  have  responses  most 
closely  matching  those  of  the  actual  aircraft." 

The  parameters  defining  the  resulting  equivalent  system  (frequency,  damping  ratio,  time  constants 
etc.)  rather  than  any  modes  of  the  high  order  system,  are  to  be  compared  with  the  specification 
requirements.  Guidance  as  to  how  the  contractor  shall  proceed  with  his  equivalent  system 
definition  is  addressed  in  the  proposed  MIL  Handbook,  Volume  II  of  reference  (f). 

The  Naval  Air  Development  Center,  as  part  of  its  effort  in  identifying  flying  qualities  criteria 
for  manned  aircraft, undertook  the  determination  of  equivalent  system  descriptions  of  current  Navy 
tactical  aircraft.  The  determination  of  classical  pitch-rate  short-period  models  for  current  fleet 
aircraft  was  reported  in  reference  (g).  This  report  presents  similar  results  for  the  lateral-directional 
responses  of  the  same  U.S.  Navy  aircraft. 

PURPOSE 

The  purpose  of  this  effort  was  to  investigate  the  utility  of  the  equivalent  systems  approach  to 
defining  the  dynamic  lateral-directional  flying  qualities  parameters  of  augmented  aircraft.  This 
report  presents  equivalent  low  order  system  models  for  current  U.S.  Navy  tactical  aircraft  and 
compares  them  with  the  modal  requirements  of  Ml  L-F-8785C,  reference  (e). 

SCOPE 

Lateral-directional  modal  responses  for  the  A-6,  A-7,  F-14,  F-18  and  S-3  aircraft  were 
analyzed  in  this  effort.  Where  applicable,  each  aircraft  was  assumed  to  have  its  Stability/Control 
Augmentation  System  (SAS/CAS)  ON.  The  flight  conditions  investigated  included  both  Power 
Approach  (PA)  and  Cruise  (CR)  configurations  as  presented  in  table  I. 


i . 
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TABLE  I.  FLIGHT  CONDITIONS 


Gross 

CG 

Aircraft 

Configuration 

Weight 

Position 

Altitude 

Airspeed 

(lb) 

(%  MAC) 

(ft) 

(M/KEAS) 

A-6 

CR 

39505 

23.6 

20000 

0.4 

0.72 

0.87 

S-3 

CR 

36320 

21.7 

15000 

0.36 

0.71 

F-14 

CR 

51015 

8.2 

15000 

0.40 

8.2 

0.715 

9.9 

0.795 

10.7 

0.91 

PA 

44030 

9.8 

0 

0.183/121 

A-7 

CR 

21890 

30.0 

15000 

0.3 

0.6 

0.9 

F-18 

CR 

29930 

25.0 

10000 

0.5 

METHOD 

Frequency  response  matching  techniques  were  utilized  to  determine  low  order  equivalent  sys¬ 
tems  describing  the  complex  aircraft  high  order  representations.  Digital  computer  programs,  pre¬ 
pared  by  the  McDonnell  Aircraft  Company,  utilized  a  direct  Rosenbrock  search  algorithm 
(reference  c)  to  match  a  Bode  plot  describing  the  high  order  pilot  input  to  aircraft  output  transfer 
function  with  an  equivalent  low  order  system.  Since  this  analysis  is  concerned  with  determining 
equivalent  lateral-directional  models,  the  roll  and  sideslip  angle  responses  to  pilot  control  inputs 
were  analyzed. 

In  order  to  use  the  matching  routines,  a  description  of  the  frequency  response  of  the  system 
to  be  matched  is  required.  This  may  be  in  the  form  of  either  (1)  transfer  functions  or  (2)  numerical 
phase-gain  data  obtained  at  various  input  frequencies.  Since  only  limited  numerical  response  data  is 
available  for  the  subject  airplanes  (and  it  is  generally  corrupted  with  instrumentation  noise  and  air 
turbulence)  the  transfer  function  input  approach  was  chosen.  Each  aircraft's  transfer  functions 
describing  the  desired  responses  were  obtained  either  directly  from  available  information  (A-7)  or 
computed  via  NADC  transfer  function  programs  from  stability  and  control  derivative  information 
(A-6,  S-3,  F-14,  F-18).  Reference  (h)  through  (m)  were  used  to  obtain  this  information  as  well  as  a 
description  of  the  respective  control  systems.  With  the  aircraft's  unaugmented  dynamics  thus 
obtained,  the  control  components  present  in  each  aircraft's  control  system  (i.e.,  actuators,  stick 
feel  system,  feedback  loops,  compensation  networks,  etc.)  were  added  to  obtain  the  high  order 
transfer  function  describing  each  aircraft/control  system  combination  and  flight  condition.  Brief 
descriptions  of  the  aircraft  and  their  respective  control  systems  are  presented  in  appendix  A. 
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The  equivalent  low  order  systems  were  obtained  via  a  frequency  response  matching  technique 
which  minimizes  a  mismatch  function.  This  mismatch  function  is  defined  as  the  weighted  sum  of 
the  squares  of  the  differences  in  magnitude  and  phase  angle  between  the  high  and  low  order  systems 
at  a  number  of  discrete  frequencies.  Quantitatively,  this  can  be  expressed  as: 


wn 

M  =  13  [(GhoS  -  Gl_OS)2  +  -01745  (4>HOS  ~  $LOS)25 

co  1 


(1) 


where  G  equals  the  gain  in  decibels  and  4>  is  the  phase  in  degrees.  Summing  the  mismatch  function 
over  a  number  of  frequencies,  evenly  spaced  on  a  logarithmic  scale,  is  similar  to  minimizing  the 
integral  of  the  square  of  the  error  on  a  Bode  plot.  As  a  result,  it  is  possible  to  qualitatively  compare 
the  matches  with  the  quantitative  mismatch  results.  The  frequency  range  over  which  the  minimiza¬ 
tion  was  conducted  was  chosen  to  span  the  pilot's  primary  frequency  range  of  interest  (0.1  —  10 
rad/sec). 


The  frequency  response  matching  procedure  enables  the  analyst  to  match  any  high  order 
system  with  any  desired  low  order  system  format.  Restricting  the  scope  of  this  analysis  to  the 
determination  of  dynamic  lateral  directional  characteristics  establishes  the  possible  forms  of  the 
order  system  to  be  matched.  Beginning  with  the  three  degree  of  freedom  equations  of  motion 
describing  the  lateral  directional  responses,  the  MIL-SPEC  dynamic  lateral  directional  characteristic:. 
(fOR-  WDR.  Tr  and  Ts>  can  be  obtained  from  the  transfer  functions  relating  roll  angle  and  sideslip 
angle  response  to  control  deflection: 


0(5)  _  Lsa  (s2  +  <t>2) _ 

^a(s)  (s+1/rr)  (s+1/ts)  (s2  +  2f[)pa;Q)ps  +  u>[}R2) 

ff(s)  N5r  (s-H/t^)  (s+1/r02)  (s+1/T0g) 

(s+1/rr)  (s+1/ts)  (s2  +  2fDRojQps  +  ojqr2) 


Equations  (2)  and  (3)  can  be  further  simplified  by  decoupling  the  lateral  and  directional  modes 
of  motion.  Assuming  only  a  single  degree  of  freedom  roll  response,  the  roll  equation  of  motion 
yields  the  following  approximate  transfer  function: 


0(s)  _  *~^a 

aa(s)  “  s  (s+1/rr) 


(4) 


Similarly  by  eliminating  the  rolling  degree  of  freedom,  the  Dutch  roll  approximation  may  be  ex¬ 
pressed  as: 


m  =  Nfir _ 

(s2  +  2?dro>drs  +  o;dR2) 


(5) 


In  matching  the  response  of  high  order  and/or  closed  loop  control  systems  to  pilot  command 
inputs,  the  control  input  required  in  equations  (2)  through  (5)  is  the  cockpit  control  force  or 
deflection. 
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The  inclusion  of  control  system  components  in  the  high  order  system  may  result  in  high 
frequency  phase  responses  which  exceed  the  phase  asymptotes  of  the  classical  systems.  Therefore, 
the  equivalent  system  procedure  augments  equations  (2)  through  (5)  via  the  addition  of  a  time 
delay.  This  delay  introduces  phase  lag  without  altering  the  gain  characteristics.  The  resulting 
equivalent  systems  may  then  be  expressed  as  follows: 

Complete  three  degree  of  freedom 


— ti  s 

<fi(s)  _  K0  (s2  +  2^u?^s  +  cj02)  e _ 

(s+1/rr)  (s+1/ts)  (s2  +  2foRGjQRS  +  ojdr2) 

(s+1/r^1 )  (s+1  /r^2>  (s+l/r^)  e 
^rp^s)  (s+1/Tr)  (s+1/ts)  (s2  +  2fQpojQRS  +  cjqR^ 

Approximate  forms: 

d>(s)  =  H  ^ 

8ap(s)  s  (s+1/rr) 

Pis)  =  « P  e~^S _ 

5rp<s)  (s2  +  2fDRWQRS  +  cjdR2) 


(6) 


(7) 


(8) 

19) 


In  order  to  further  simplify  the  approximate  forms,  roll  rate,  rather- than  roll  angle,  was  used 
to  evaluate  the  roll  response.  The  approximate  roll  rate  response  appears  as: 


Si s)  = 

5ap(s)  (s+1/rr) 


(10) 


Equations  (6)  and  (7)  have  been  implemented  by  Hodgkinson  in  a  digital  computer  program 
called  LATFIT,  reference  (n).  This  program  allows  for  either  individual  or  simultaneous  matching 
of  the  roll  angle  and  sideslip  responses.  When  using  the  LATFIT  program,  initial  estimates  for  the 
unknown  parameters  must  be  supplied  by  the  analyst. 


It  should  be  noted  that  the  guidance  for  determining  equivalent  system  models  given  in  refer¬ 
ence  (f)  is  to  individually  match  the  approximate  sideslip  response  via  equation  (9)  to  obtain  Dutch 
roll  information  and  the  complete  roll  rate  response  via 


p(s) 


=  s  =  s-  equation  (6) 


(ID 
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to  obtain  roll  mode  information.  It  will  be  shown  in  the  Results  and  Discussion  section  that  it  is 
often  necessary  to  match  the  sideslip  and  bank  angle  responses  simultaneously,  in  which  case 
consistent  models  should  be  utilized. 

In  per.  jrming  the  equivalent  system  match,  the  analyst  worked  interactively  with  the  program 
to  determine  which  of  the  decision  variables  to  vary  in  the  search  for  the  equivalent  low  order  sys¬ 
tem.  In  general,  the  procedure  outlined  below  was  utilized: 

a)  Initially,  the  0  and  j3  response  transfer  functions  were  matched  independently  with  the 
approximate  forms  of  the  low  order  equivalent  system.  (This  step  was,  in  general,  easily 
accomplished  and  provided  qualitatively  acceptable  matches  when  compared  in  both  the 
frequency  and  time  domains.) 

b)  Secondly,  improvement  in  the  mismatch  function  was  evaluated  by  obtaining  the  com¬ 
plete  low  order  forms,  equations  (6)  and  (7).  This  was  accomplished  by  starting  foR, 
wdr,  rr,  and  the  time  delays  from  the  approximate  form  results.  In  addition 

and  T0 2  were  either  started  at  their  high  order  system  values  or  at  values  close  to 
?DR,  wdr,  and  rr,  respectively.  The  remaining  factors,  ,  t ^  and  Ts.  which  were 

generally  outside  the  frequency  range  of  interest,  were  held  fixed  at  their  high  order 
system  values. 

The  proximity  of  the  roll  angle  numerator  and  denominator  oscillatory  roots,  as  well  as  one 
of  the  sideslip  numerator  time  constants  (r^)  and  the  roll  mode  time  constant,  restricted  the  low 
order  equivalent  system  determination  when  matching  0  and  0  independently  with  the  complete 
form  (step  b).  This  problem  could  be  alleviated  by  fixing  the  numerator  roots  at  their  high  order 
system  value  or  by  simultaneously  matching  0  and  0  while  constraining  the  equivalent  denomina¬ 
tors  to  be  identical.  While  the  former  case  generally  produced  satisfactory  results,  knowledge  of 
the  numerator  roots  may  not  always  be  available  prior  to  the  matching  process.  In  the  latter  case 
the  roll  mode  time  constant,  identifiable  from  the  roll  response,  helped  to  locate  the  sideslip 
numerator  root  while  the  Dutch  roll  roots,  identifiable  from  the  sideslip  response,  helped  to 
locate  the  roll  angle  numerator  roots. 

When  individual  control  system  roots  were  present  in  the  frequency  range  of  interest,  as  with 
prefilters  for  example,  the  problems  associated  with  determining  an  acceptable  equivalent  system 
match  were  compounded. 

The  following  technique  (further  explained  in  the  Results  and  Discussion  section)  was  found 
to  overcome  these  problems  with  a  minimum  amount  of  computational  effort: 

a)  Determine  an  approximate  denominator  root  location  using  equations  (9)  and  (10). 

b)  Simultaneously  match  equations  (6)  and  (7)  in  the  following  sequence  (Fix  r r^,  and 
rs  at  their  corresponding  high  order  system  values  in  these  runs): 

1.  Fix  Cdr,  cjqr,  and  Tr  at  the  values  obtained  in  step  a.  Match  all  remaining 
parameters  (beginning  with  «  ?qR'  w0  ^  ^DR-  and  Tfc  55  T r>- 

2.  Fix  $0,  cu0,  and  0^  at  the  values  obtained  in  step  b.1.  Start  all  other  parameters  from 
their  step  b.  1  values  and  obtain  a  new  match. 
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At  this  point,  all  the  parameters  could  be  started  at  the  results  of  step  b.2  and  freed  in  the  matching 
algorithm  to  obtain  a  final  match.  However,  the  matching  technique  did  not  always  provide  a 
unique  solution.  In  addition,  the  computational  effort  expended  was  not  justified  by  the  small 
mismatch  improvement  generally  obtained.  Therefore,  the  technique  utilized  was  to  terminate 
the  matching  procedure  after  step  b.2  and  compare  those  res  'ting  equivalent  parameters  with  the 
MIL-F-8785C  requirements. 

In  performing  these  matches,  the  "goodness  of  fit"  was  determined  quantitatively  by  the  value 
of  the  mismatch  parameter  and  qualitatively  from  Bode  plots  and  time  history  responses  to  unit 
impulse  and  step  control  inputs.  Although  no  minimum  value  of  mismatch  has  been  conclusively 
established,  several  guidelines  have  been  proposed.  Hodgkinson  (reference  c),  initially  recom¬ 
mended  that  a  good  match  is  obtained  when  the  mismatch  function  is  less  than  or  equal  to  10.  In 
a  limited  experimental  effort  by  Smith  (reference  o),  pilot  ratings  were  correlated  with  various 
levels  of  mismatch.  These  results  indicated  that  pilots  were  insensitive  to  analytical  mismatches  up 
to  200  in  the  frequency  range  of  0.1  to  10.0  radians/second. 
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RESULTS  AND  DISCUSSION 


GENERAL 

Approximate  and  complete  lateral  directional  equivalent  system  models  were  determined  for 
the  A-6,  A-7,  S-3,  F-14,  and  F-18  airplanes.  Each  of  these  airplanes  incorporates  rate  and/or 
acceleration  feedback  for  stability  augmentation.  The  A-7,  F-14  and  F-18  airplanes  include  for¬ 
ward  loop  control  augmentation.  Similar  techniques  were  utilized  to  obtain  the  equivalent  system 
matches  for  each  of  these  airplanes. 

In  general,  acceptable  equivalent  system  models  were  developed  for  each  of  these  aircraft.  The 
greatest  difficulties  were  experienced  for  those  configurations  with  forward  loop  prefilter  like 
components  whose  frequencies  lie  in  the  range  of  0.1  -  10  rad/sec.  The  purpose  of  these  prefilters 
is  to  provide  attenuation  of  high  frequency  inputs.  The  overall  effect  is  to  impose  both  an  attenua¬ 
tion  and  a  lag  on  the  system  response.  The  equivalent  system  model  reflects  this  effect  via  a  modi¬ 
fied  modal  response  parameter  and  a  significant  increase  in  time  delay. 

Comparisons  of  the  equivalent  system  parameters  with  the  specification  requirements  of  refer¬ 
ence  (c)  resulted  in  Level  1  flying  qualities  for  all  but  a  few  cases.  At  .3M,  the  A-7  equivalent  roll 
mode  time  constant  exceeds  the  Level  1  boundary.  Also,  the  F-18  configuration  investigated  ex¬ 
ceeds  the  roll  numerator  time  delay  requirement  of  0.10  seconds. 

EQUIVALENT  SYSTEM  MODELS 

The  emphasis  in  lateral  directional  equivalent  systems  matching  has  been  centered  on  the  roll 
and  sideslip  responses,  the  parameters  of  which  are  necessary  for  comparison  with  the  require 
ments  of  Ml  L-F-8785C.  Initial  matches  were  obtained  for  the  roll  and  sideslip  angle  responses  to 
pilot  control  input  transfer  functions  for  the  S-3  and  A-6.  In  general,  these  configurations  were 
characterized  by  the  bare  airframe  augmented  by  control  surface  actuators  and  rate  feedbacks. 

From  these  configurations,  a  basic  understanding  of  the  frequency  response  matching  process  was 
obtained.  Subsequently,  the  effect  of  compounding  dynamic  components  such  as  lateral  accelera¬ 
tion  feedback  and  forward  loop  lags  were  investigated  for  the  A-7,  F-14,  and  F-18  aircraft. 

S-3  Airplane 

The  S-3  airplane's  lateral  directional  control  system  utilizes  feedback  from  the  yaw  rate  gyro 
to  augment  the  aircraft's  basic  Dutch  roll  characteristics.  The  resulting  roll  and  sideslip  angle 
responses  to  pilot  control  inputs  can  be  represented  by  fifth  and  sixth  order  numerators,  respec¬ 
tively,  over  an  eighth  order  denominator.  The  feedback  components  add  a  numerator  and  denomi¬ 
nator  root  in  the  vicinity  of  0.35  rad/sec  while  the  actuator  adds  a  single  denominator  root  at  ap¬ 
proximately  22  rad/sec  in  the  lateral  axis  and  29  rad/sec  in  the  directional  axis.  Because  of  the 
proximity  of  the  added  feedback  roots,  they  effectively  cancel  each  other  while  the  actuator  roots 
are  outside  the  pilot's  frequency  range  of  interest.  Approximate  pole-zero  cancellation  of  roots 
within  the  frequency  range  of  interest  results  in  a  transfer  function  similar  in  form  to  that  of  the 
desired  equivalent  system.  It  is  therefore  expected  that  a  close  match  of  the  higher  order  system 
will  be  obtained. 

The  McAir  recommended  matching  procedure  was  utilized  with  the  NAVFIT  and  LATFIT 
computer  programs  to  determine  an  equivalent  system  for  each  of  the  flight  conditions  analyzed. 
Both  the  approximate  and  complete  forms  of  the  roll  and  sideslip  angle  transfer  functions  were 
used.  The  following  discussion  will  focus  on  a  representative  case  analyzed  for  the  S-3  aircraft  to 
demonstrate  the  matching  procedures. 
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The  equivalent  system  parameters,  along  with  the  high  order  system  values,  for  the  S-3  air¬ 
craft  at  15,000  ft  and  0.36M  are  shown  in  table  II.  The  approximate  forms  (trials  1  and  2)  yield 
good  matches  of  roll  mode  time  constant  and  Dutch  roll  frequency  and  damping  ratio  when  com¬ 
pared  to  the  high  order  system  values.  Equivalent  time  delays  are  acceptable  and  mismatch  values 
are  relatively  low.  The  frequency  and  time  responses  for  the  high-order  and  low-order  equivalent 
systems  are  compared  in  figures  1  and  2.  As  indicated  in  these  figures,  the  frequency  responses  are 
relatively  well  matched  in  the  region  from  .1  to  10  rad/sec.  This  impression  is  confirmed  by  the 
excellent  match  of  rise  time  on  roll  rate  (roll  mode  time  constant  and  control  power  combination) 
and  the  oscillatory  response  in  sideslip  angle  (frequency  and  damping).  There  are  however,  minor 
differences  between  the  high-and  low-order  responses;  the  equivalent  steady  state  response  does  not 
decay  with  time  (a  result  of  the  large  mismatch  exhibited  at  frequencies  less  than  0.1  rad/sec)  and 
the  oscillatory  roll  response  is  absent  in  the  low-order  system.  These  differences  result  from  their 
not  being  modelled  in  the  approximate  forms.  Matching  the  complete  forms  should  improve  both 
of  these  areas. 


Improvements  in  the  steady  state  response  can  be  demonstrated  by  adding  the  low  frequency 
roots  to  the  equivalent  system  models  and  identifying  the  roots  over  an  expanded  frequency  range. 
The  equivalent  models  can  then  be  represented  as: 


0  _  K<t>e  ^ 

(sVl/rr)  (s  +  1/rs) 

0  K0IS+  1/r^)  e-t0S 

^PED  (s  +  1/rs)  (s^  +  2  codr  s  + 


(ID 


(12) 


Expanding  the  frequency  match  range  to  encompass  the  added  low  frequency  breakpoints,  the 
results  presented  in  trials  3  and  4  of  table  II  were  obtained.  In  the  case  of  the  roll  response,  the 
NAVFIT  program  was  used,  with  arbitrary  initial  parameter  values,  to  obtain  the  results  shown. 
However,  in  the  case  of  the  sideslip  response,  it  was  necessary  to  utilize  the  identified  roll  mode 
spiral  root  as  a  starting  point,  in  the  NAVFIT  program,  in  order  to  obtain  satisfactory  results. 

The  mismatch  between  the  high-and  low-order  system  roll  response  was  greatly  improved, 
as  shown  in  figure  3.  The  only  remaining  discrepancy  between  the  two  systems  is  now  the  absence 
of  the  oscillatory  component  in  the  low  order  response.  The  differences  in  the  high-  and  low-order 
sideslip  response  are  now  almost  inperceptible,  as  shown  in  figure  4. 

Having  demonstrated  the  ability  to  identify  the  low  frequency  roots,  and  the  associated  mis¬ 
match  improvement,  rs  and  t p.  were  fixed  at  their  known  high-order  values  in  all  subsequent  runs 
in  order  to  simplify  the  matching  process. 

The  complete  form  equations  were  used  in  an  attempt  to  further  improve  the  response 
matches  in  the  mid  frequency  range.  In  this  case,  independently  matching  the  roll  or  sideslip  angle 
responses  was  difficult  due  to  the  proximity  of  several  numerator  and  denominator  roots.  The 
oscillatory  roll  angle  response  roots  [?0,  W0]  and  t fD R*  ^DRl  were  close  enough  to  prevent  the 
search  algorithm  from  identifying  realistic  values.  Similar  difficulties  were  also  encountered  for 
the  sideslip  angle  response  due  to  the  proximity  of  a  sideslip  numerator  constant  (r^)  and  the  roll 
mode  time  constant  (rr).  This  difficulty  was  overcome  (trials  5  and  6  of  table  II)  by  fixing  the 
parameters  f0,  cj0,  and  rp^  at  their  high  order  system  values  or  starting  them  very  close  to  these 
values.  The  resulting  match  statistics  indicate  a  decrease  in  both  the  roll  and  sideslip  mismatch. 
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(Note:  See  Appendix  B  for  summary  of  all  flight  conditions.) 
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ure  3.  S-3  Roll  Response  -  Approximate  Form  Plus  Low  Frequency  Root  (Trial  3) 
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Also  evident  is  that  the  roll  mode  time  constant  identified  for  each  case  closely  agrees  with  the  high 
order  system  value.  This  technique,  however,  requires  apriori  knowledge  of  the  high  order  system 
modal  parameters,  which  may  not  always  be  available. 

The  lack  of  apriori  knowledge  of  the  high  order  system  numerator  roots  can  be  avoided  by 
either  1)  fixing  the  denominator  values  at  the  identified  approximate  form  results  and  performing 
the  required  match,  or  2)  simultaneously  matching  roll  and  sideslip  angle  responses  while  constrain¬ 
ing  the  equivalent  denominators  to  be  identical.  The  first  procedure  is  demonstrated  in  trials  7 
and  8  of  table  II.  The  results  obtained  show  only  slight  variations  from  those  obtained  by  identify 
ing  denominator  roots  with  the  numerator  roots  held  fixed.  Then  beginning  with  the  roots  at  these 
latest  locations,  it  was  possible  to  free  all  parameters,  as  demonstrated  in  trials  9  and  10  of  table  II. 

In  the  second  procedure,  the  Dutch  roll  roots,  easily  identified  from  the  sideslip  response, 
helped  to  locate  the  roll  angle  numerator  roots  while  the  roll  mode  time  constant,  easily  identified 
from  the  roll  angle  response,  helped  to  locate  the  sideslip  numerator  root.  The  results  for  this 
procedure  are  presented  in  the  final  column  of  table  II.  Excellent  agreement  was  obtained  between 
all  high-  and  low-order  modal  parameters  investigated  along  with  extremely  low  values  of  mismatch. 
Frequency  and  time  responses  for  the  simultaneous  matching  procedure  are  presented  in  figures  5 
and  6.  Excellent  agreement  is  shown  between  the  high-  and  low-order  systems  with  the  oscillatory 
component  of  the  roll  response  now  being  identified. 

The  approximate  and  complete  form  sideslip  angle  results  provide  additional  insight  into  the 
time  delay  parameter.  In  the  approximate  form,  the  time  delay  accounts  for  the  total  phase  lag 
imparted  by  high  frequency  roots,  whereas  in  the  complete  system,  a  high  frequency  aircraft 
numerator  root  ( l/r^)  is  included  in  the  equivalent  model.  Since  this  numerator  root  adds  high 
frequency  phase  lead,  the  complete  form  time  delay  has  been  increased  in  order  to  match  the  over 
all  high  frequency  phase  characteristics.  The  approximate  time  delay  can  therefore  be  interpreted  as 
representing  the  overall  aircraft/control  system  lag  experienced  by  the  pilot,  while  the  complete 
system  time  delay  represents  the  time  delay  contribution  arising  from  the  control  system. 

The  equivalent  system  models  obtained  for  the  S-3  airplane  flight  conditions  investigated  are 
summarized  in  Appendix  B. 

A-6  Airplane 

The  A-6  airplane  utilizes  both  roll  rate  and  yaw  rate  feedbacks  to  augment  the  aircraft's  basic 
lateral -directional  stability  characteristics.  The  resulting  control  system  descriptions  add  numerator 
and  denominator  roots  in  the  vicinity  of  0.2  to  0.5  rad/sec  while  the  flaperon  and  rudder  actuators 
add  denominator  roots  at  approximately  19.6  and  27.0  rad/sec,  respectively. 

The  following  discussion  is  for  a  representative  flight  condition  (0.40M,  20,000  ft  altitude)  for 
the  A-6  airplane.  Results  for  all  flight  conditions  analyzed  are  summarized  in  Appendix  B. 

Transfer  functions  describing  the  airplane's  roll  and  sideslip  angle  responses  to  pilot  control 
inputs  can  be  represented  by  fifth  and  sixth  order  numerators,  respectively,  over  an  eighth  order 
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denominator.  Inspection  of  the  high-order  system  roots  in  the  frequency  range  of  interest  yields 
the  following  observations: 

a)  Utilizing  approximate  pole-zero  cancellation,  only  one  significant  numerator,  and  two 
significant  denominator  roots  remain  in  the  roll  rate  transfer  function.  The  effect  of 
these  roots  is  to  provide  magnitude  attenuation  of  20  db  per  decade  of  frequency  and  a 
high  frequency  phase  lag  approximately  90  deg  less  than  that  at  low  frequencies.  Since 
these  are  the  same  response  characteristics  exhibited  by  the  low  order  system  model  it 
can  be  expected  that  an  equivalent  system  match  can  !’e  readily  determined. 

b)  Applying  approximate  pole-zero  cancellation  to  the  sideslip  response,  only  a  single  pair 
of  oscillatory  denominator  roots  remains  in  the  frequency  range  of  interest.  Since  the 
form  of  the  low  order  system  is  identical  to  this  representation  of  the  high  order  system, 
an  excellent  sideslip  angle  match  is  anticipated. 

These  observations  are  borne  out  in  the  equivalent  system  results  presented  in  table  III. 

Good  matches  of  both  the  roll  rate  and  sideslip  angle  responses  via  the  approximate  forms 
were  obtained,  as  shown  in  figures  7  and  8. 

The  complete  form  matching  technique  results  are  similar  to  those  obtained  for  the  S-3  air¬ 
plane.  In  the  case  of  the  independent  response  parameter  matching  technique  it  was  again 
necessary  to  fix  fa  and  at  their  high  order  system  values.  Significant  differences  are  noted 
between  the  bank  and  sideslip  angle  results  for  rr  and  {"qr,  while  cuqr  shows  only  minor  differences. 

In  order  to  obtain  a  consistent  set  of  modal  parameters,  it  was  necessary  to  simultaneously 
match  both  bank  and  sideslip  angle  responses,  while  constraining  the  identified  denominators  to  be 
identical.  The  numerator  parameters,  fa,  W0  and  r ^  were  again  allowed  to  go  free  in  the  matching 
process.  The  results,  presented  in  table  III,  are  seen  to  be  a  compromise  on  the  individually  deter¬ 
mined  low-order  systems  —  the  mismatch  values  are  increased  and  the  denominator  parameters, 

Tr,  {“dr,  and  u>dr  lie  near  the  values  obtained  for  the  independent  matching  technique.  Figures 
9  and  10  show  the  frequency  and  time  responses  for  these  simultaneous  matches. 

The  mistmatch  values  for  both  the  approximate  and  complete  forms  are  nigher  than  those  ob¬ 
tained  for  the  S-3  airplane.  The  question  of  acceptable  mismatch  was  addressed  by  Smith, 
Hodgkinson,  and  Snyder  (reference  p).  Based  on  limited  experimental  data,  they  concluded  that 
pilot  ratings  were  insensitive  to  analytical  mismatch  values  up  to  190.  Since  the  maximum  value 
obtained  for  the  A -6  airplane  was  145.4,  these  matches  were  considered  to  be  acceptable.  Figures 
7  and  8  show  the  frequency  and  time  responses  for  the  approximate  forms  and  graphically  illus¬ 
trate  these  levels  of  mismatch. 

A-7  Airplane 

The  A-7  airplane  utilizes  forward  path  compensation  as  well  as  yaw  and  roll  rate  feedback,  as 
described  in  Appendix  A,  to  augment  the  aircraft's  basic  stability  characteristics.  The  roll  and  side¬ 
slip  angle  responses  to  pilot  control  inputs  can  be  represented  by  seventh  and  tenth  order  numera¬ 
tors,  respectively,  over  twelfth  order  denominator  transfer  functions.  The  forward  mechanical  path 
includes  feel  system  components  which  introduce  first  order  roots  at  approximately  12  rad/sec. 

The  electrical,  or  command  augmentation,  feed  forward»path  includes  prefilter  type  components  at 
3  and  10  rad/sec.  The  3  rad/sec  root  lies  well  within  the  frequency  range  of  interest  and  could  be 
expected  to  cause  problems  in  obtaining  an  equivalent  system  match.  The  feedback  components 
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TABLE  ill 

A-6  AIRPLANE  EQUIVALENT  TRANSFER  FUNCTIONS 
Roll  and  Sideslip  Angle  Response  to  Cockpit  Control  Inputs 
CR  20,000  FT  .40  M 


HOS 

VALUES 

APPROXIMATE  FORM 

COMPLETE  FORM 

INDEPENDENT 

INDEPENDENT 

SIMULTANEOUS 

P 

P 

0 

P 

<p&P 

17.84 

.670 

- 

.702 

- 

.646 

.283 

- 

- 

.283* 

- 

.194 

w0 

1.77 

- 

- 

1.77* 

- 

1.56 

<0 

- 

.0076 

- 

.018 

- 

.0077 

K/3 

.0078 

- 

.0293 

- 

.0003 

.0003 

% 

131.58 

- 

- 

- 

-131.58* 

-131.58* 

r02 

.148 

— 

— 

— 

.172 

.692 

T03 

.010 

— 

— 

— 

.010* 

.010* 

XP 

- 

- 

.025 

- 

.034 

.042 

rr 

.156 

.461 

- 

.767 

.181 

.642 

Ts 

105.26 

- 

- 

105.26* 

105.26* 

105.26* 

?DR 

.269 

- 

.251 

.465 

.258 

.299 

cc 

Q 

3 

1.73 

- 

1.736 

2.17 

1.74 

1.71 

M , 

0 

- 

145.4 

- 

101.0 

- 

121.8 

- 

4.4 

- 

.46 

4.2 

'Fixed  parameter 


(Note:  See  Appendix  B  for  summary  of  all  flight  conditions) 
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add  numerator  and  denominator  roots  in  the  vicinity  of  both  1  and  31  rad/sec,  which  approxi¬ 
mately  cancel  each  other.  In  addition,  the  surface  actuators  add  denominator  roots  at  approxi¬ 
mately  20  rad/sec,  which  are  outside  the  frequency  range  of  interest. 

The  following  discussion  is  for  a  representative  flight  condition  at  0.6  Mach  and  15,000  ft 
altitude.  The  results  for  all  A-7  conditions  analyzed  are  summarized  in  Appendix  B.  This  dis¬ 
cussion  is  divided  into  the  results  with  and  without  the  control  augmentation  path  prefilters,  as 
significant  differences  were  obtained. 

Table  IV  summarizes  the  results  obtained  for  the  A-7  aircraft  at  15,000  ft  and  0.60  M  without 
the  control  augmentation  path  prefilters.  Using  the  approximate  forms  for  the  low  order  system, 
good  match  statistics  were  obtained  when  compared  to  the  high  order  system  values.  Roll  mode 
time  constant  shows  excellent  agreement,  while  Dutch  roll  frequency  and  damping  ra/o  show  only 
minor  differences.  The  mismatch  function  was  relatively  low  for  both  the  roll  and  sideslip  angle 
responses. 

Frequency  and  time  response  comparisons  for  the  approximate  forms  are  presented  in  figures 
1 1  and  12.  The  only  significant  differences  occur  in  the  steady  state  portions  of  the  responses 
which  may  be  attributed  to  the  lack  of  a  spiral  mode  root  in  the  equivalent  systems. 

The  complete  form  results  are  also  shown  in  table  IV.  Again,  and  had  to  be  either 

fixed  at  or  started  very  close  to  their  high  order  system  values.  If  these  parameters  were  not  fixed, 
their  resulting  values  were  not  unique  but  rather  varied  as  their  initial  value  varied.  This  was 
especially  true  of  the  roll  angle  response  as  show  in  table  V.  Unique  values  for  r ^  and  rr  could  be 

determined  when  matching  sideslip  angle  along,  but  only  after  a  very  large  number  of  iterations 
had  been  completed. 

This  difficulty  was  overcome  by  simultaneously  matching  the  roll  and  sideslip  angle  responses 
as  was  previously  done  for  the  S3  and  A-6  aircraft.  Table  IV  presents  the  results  of  this  simul¬ 
taneous  match  for  the  A-7  and  indicates  good  match  statistics.  The  mismatch  functions  are  very 
low  while  the  modal  parameters  are  very  close  to  the  high  order  system  values.  Frequency  and  time 
history  responses  for  these  results  are  shown  in  figures  13  and  14.  There  is  essentially  no  difference 
between  the  high  order  and  the  low  order  equivalent  systems.  The  steady  state  responses  identified 
with  the  simultaneous  matching  technique  show  a  significant  improvement  when  compared  to  the 
approximate  form  (figures  11  and  12). 

The  inclusion  of  the  roll  command  augmentation  in  the  forward  loop  had  a  significant  impact 
on  both  the  commanded  response  and  the  equivalent  system  model.  The  purpose  of  the  prefilters  is 
to  provide  attenuation  of  inputs  at  frequencies  greater  than  approximately  3  rad/sec.  The  effect  of 
this  attenuation  is  to  impose  a  lag  in  the  roll  response  to  step  control  inputs  as  shown  in  figure  15. 

Approximate  pole-zero  cancellation  of  the  high  order  bank  angle  response,  in  the  frequency 
range  of  interest,  results  in  a  second  over  fifth  order  transfer  function.  Since  this  representation  is 
different  from  that  of  the  equivalent  models,  difficulties  in  obtaining  an  equivalent  system  match 
can  be  anticipated.  The  additional  phase  lag  introduced  by  the  command  augmentation  prefilters 
is  accounted  for  in  the  equivalent  system  model  via  the  time  delay.  However,  since  there  are  also 
gain  modifications  resulting  from  the  prefilters  in  the  frequency  range  of  interest,  which  are  not 
accounted  for  in  the  equivalent  system  model,  the  resulting  equivalent  system  parameters  are,  at 
best,  an  average  characteristic  reflecting  the  contribution  of  aircraft  and  control  system  com¬ 
ponents.  This  condition  is  reflected  in  the  equivalent  system  parameters  for  the  A-7  aircraft  as 
shown  in  table  VI. 
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TABLE  IV 

A-7  AIRPLANT  EQUIVALENT  TRANSFER  FUNCTIONS,  NO  PREFILTERS 
Roll  and  Sideslip  Angle  Response  to  Cockpit  Control  Inputs 
CR  15,000  FT.  0.6  M 
CONTROL  AUGMENTATION  OFF 


HOS 

VALUES 

APPROXIMATE 

FORM 

COMPLETE 

FORM 

INDEPENDENT 

IND 

EPENDENT 

SIMULTANEOUS 

0 

(3 

0 

0  &  & 

K0 

422.2 

23.2 

- 

21.25 

- 

20.16 

b 

.56 

- 

- 

.56* 

- 

.449 

^0 

2.02 

- 

- 

2.02* 

1.92 

.053 

- 

.045 

.040 

.0011 

- 

.0061 

- 

.000066 

.000065 

r0i 

-1250.0 

- 

- 

- 

-1250.0* 

-1250.0* 

tP2 

.198 

- 

- 

- 

.195 

.408 

t03 

.009 

- 

- 

- 

.009* 

.009* 

— 

- 

.032 

- 

.060 

.058 

m 

.20 

.20 

- 

.23 

.14 

.25 

m 

37.45 

- 

- 

37.45* 

37.45* 

37.45* 

?DR 

.49 

- 

.40 

.49 

.38 

.39 

WDR 

2.03 

- 

2.29 

2.06 

.  . 

2.15 

2.01 

. 

- 

20.1 

- 

0.35 

- 

0.86 

Mp 

- 

10.4 

7.0 

4.1 

'Fixed  parameter 

(Note:  See  Appendix  B  for  summary  of  all  flight  conditions) 
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TABLE  V 

EFFECT  OF  INITIAL  PARAMETER  VALUES  ON  ROLL  ANGLE  MATCHING 
A-7  Airplane  0.6  Mach,  15,000  FT  Altitude 
Control  Augmentation  Off 


HOS 

APPROXIMATE 

COMPLETE 

INITIAL 

FINAL 

INITIAL 

FINAL 

K  0 

422.2 

23.2 

5.0 

19.0 

5.0 

17.4 

4 

.56 

- 

.20 

1.26 

.10 

1.60 

2.02 

- 

3.0 

2.67 

1.0 

4.50 

■ 

.20 

.20 

.25 

.59 

.25 

.38 

?DR 

.49 

- 

.35 

.96 

.35 

1.22 

WDR 

2.03 

- 

2.5 

4.20 

2.5 

5.39 

V 

- 

.053 

.05 

.036 

.05 

.032 

M0 

- 

20.1 

- 

0.39 

- 

0.35 
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Figure  15.  Impact  of  Roll  Command  Augmentation  on  A-7  Roll  Rate  Response 

Referring  first  to  the  approximate  form  results,  it  can  be  seen  that  the  time  delay  and  roll 
mode  time  constant  are  considerably  larger  than  those  for  the  Command  Augmentation  OFF  results 
of  table  IV.  Likewise,  the  bank  angle  mismatch  has  increased  significantly,  reflecting  the  difficulty 
in  obtaining  a  good  match  with  the  additional  prefilter  roots  present  in  the  high  order  system.  The 
sideslip  response,  which  is  unaffected  by  the  lateral  command  augmentation,  yields  the  same  equiva- 
lent  system  as  shown  in  table  IV. 

The  complete  form  bank  angle  results  (trial  3  of  table  VI)  again  illustrate  the  difficulties  ex¬ 
perienced  in  obtaining  equivalent  system  results  with  a  closeness  of  numerator  and  denominator 
roots.  Although  very  excellent  match  statistics  are  obtained  (M0  =  2.7),  the  modal  parameters 
identified  are  completely  unrealistic.  It  was  again  necessary  to  either  fix  the  numerator  roots  at 
their  known  higher  order  system  value  (Trial  4),  or  iteratively  fix  first  the  denominator  roots  at  the 
approximate  form  results  (Trial  5)  and  then  fix  the  resulting  numerator  roots  (Trial  6).  This  latter 
method  was  found  to  be  1 ),  the  most  straightforward  approach,  assuming  that  no  information  con 
cerning  the  root  location  was  available  apriori  and  2),  give  the  most  consistent  results.  Independent 
matching  of  the  sideslip  response  was  again  easily  accomplished,  although  the  identified  roll  mode 
time  constant  was  inconsistent  with  any  of  the  other  results.  It  should  be  possible  to  improve  the 
roll  mode  time  constant  identified  in  the  sideslip  response  by  again  iteratively  fixing  and  freeing  the 
denominator  and  numerator  terms  respectively.  However,  since  it  has  already  been  shown,  for  the 
S-3  and  A -6  airplanes,  that  simultaneous  matching  of  bank  and  sideslip  angles  simplifies  the  match¬ 
ing  procedure,  no  additional  independent  matches  were  performed. 

Freeing  all  of  the  modal  parameters  in  the  simultaneous  technique  again  provided  the  lowest 
mismatch  statistics  but  resulted  in  unrealistic  modal  parameters  (Trial  8  of  table  VI).  It  was  again 
necessary  to  iteratively  fix  and  free  the  denominator  (beginning  with  the  approximate  form  results) 
and  numerator  terms,  respectively,  in  order  to  obtain  the  equivalent  system  model  (Trials  9  and  10). 

Frequency  and  time  history  comparisons  for  the  high  order  and  equivalent  system  models, 
with  Command  Augmentation  ON,  are  presented  in  figures  16  and  17.  The  differences  evident  in 
these  two  figures  reflect  the  impact  of  not  explicitly  including  the  control  system  roots  in  the 
equivalent  system  modelling  process.  It  would  be  possible  to  add  the  additional  roots  in  the  equiva¬ 
lent  system  model  (see  for  example,  reference  (g)).  However,  the  resulting  modal  parameters  would 
not  be  consistent  with  the  comparison  data  base  used  to  generate  the  specification  requirements  of 
references  (a)  and  (e) . 
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F-14  Airplane 

The  lateral  directional  control  system  for  the  F-14  aircraft  incorporates  roll  rate,  yaw  rate  and 
lateral  acceleration  feedbacks  as  well  as  a  lateral  stick  to  differential  stabilizer  feed  forward  com¬ 
mand  to  augment  the  airplane's  basic  response  characteristics.  As  shown  in  Appendix  A,  the  lateral 
augmentation  patch  incorporates  a  prefilter  at  a  frequency  of  2.0  rad/sec.  Initially,  this  prefilter 
may  be  expected  to  impact  the  determination  of  an  equivalent  model  in  much  the  same  way  that 
the  A-7's  command  prefilter  did.  However,  the  architecture  of  the  two  systems  is  somewhat  dif¬ 
ferent:  in  the  A-7  airplane,  the  augmented  command  signal  is  fed  to  both  the  differential  tail  and 
spoilers  while  for  the  F-14,  the  augmented  command  is  directed  only  to  the  differential  tail, 
thereby  minimizing  the  prefilter  impact  on  the  total  response  model.  The  transfer  functions  repre¬ 
senting  the  roll  and  sideslip  angle  response  to  pilot  commanded  inputs  for  the  F-14  are  seventh  and 
eighth  order  numerators,  respectively,  over  a  tenth  order  denominator. 

Cancelling  numerator  and  denominator  roots  of  similar  magnitudes  and  eliminating  those  out¬ 
side  the  0.1  to  10.0  rad/sec  frequency  range,  the  roll  transfer  function  reduces  to  a  second  over 
fourth  order  response.  Since  the  low  order  equivalent  form  is  also  a  second  over  fourth  order 
response,  good  matches  were  anticipated.  In  fact,  low  order  equivalents  of  these  transfer  functions 
were  readily  obtained.  The  equivalent  system  parameters  for  a  representative  case  of  15,000  ft  and 
0.40  Mach  are  presented  in  table  VII.  The  best  matches  were  again  obtained  when  the  simul¬ 
taneous  matching  technique  was  utilized.  In  this  instance,  it  was  not  necessary  to  fix  any  of  the 
equivalent  system  parameters  in  the  range  of  .1  to  10  rad/sec.  The  corresponding  frequency  and 
time  history  responses,  shown  in  figures  18  and  19  indicate  virtually  no  difference  between  the  high 
and  low  order  systems.  The  results  for  all  F-14  airplane  flight  conditions  analyzed  are  presented  in 
Appendix  B. 
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TABLE  VII 

F-14  AIRPLANE  EQUIVALENT  TRANSFER  FUNCTIONS 
Roll  and  Sideslip  Angle  Response  to  Cockpit  Control  Inputs 
CR  15,000  FT  .40  M 


HOS 

BUB 

COMPLETE  FORM 

VALUES 

INDEPENDENT 

INDEPENDENT 

wmihmxm 

P 

0 

<t> 

0 

4>&0 

Q 

13.19 

.683 

- 

.66 

.67 

- 

.64 

4 

.70 

- 

- 

.70* 

.59 

- 

'.73 

1.28 

- 

- 

1.28* 

1.68 

- 

1.04 

V> 

- 

.054 

- 

.048 

.049 

- 

.045 

.111 

- 

.267 

- 

.0065 

.0062 

.0062 

T0\ 

-34.48 

- 

- 

- 

- 

34.48* 

-34.48* 

r02 

.388 

- 

- 

B 

- 

.388* 

_ 

1.53 

1.935 

t03 

.02 

- 

- 

- 

- 

.02* 

.02* 

'0 

- 

- 

.020 

- 

- 

.06 

.056 

.054 

m 

.36 

.671 

- 

.55 

.53 

.23 

.55* 

.701 

m 

-62.50 

- 

- 

-62.50* 

-62.50* 

-62.50* 

^DR 

.61 

- 

.491 

.55 

.49* 

.43 

.53 

.591 

WDR 

1.07 

- 

1.515 

1.18 

1.515* 

1.31 

1.06 

1.06 

- 

12.6 

- 

1.8 

4.1 

- 

- 

1.4 

Mp 

- 

- 

38.0 

B 

- 

18.3 

2.6 

1.5 

‘Fixed  parameter 

(Note:  See  Appendix  B  for  summary  of  all  flight  conditions) 
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F-18  Airplane 

The  F-18  lateral  directional  control  system  is  part  of  a  highly  complex  digital  flight  control 
system.  Lateral  acceleration,  roll  rate  and  yaw  rate  signals  are  shaped  and  gain  scheduled  before 
being  summed  with  the  command  input  signals.  The  transfer  functions  representing  the  F-18  air¬ 
plane's  roll  and  sideslip  angle  responses  are  composed  of  20th  and  22nd  order  numerators, 
respectively,  over  a  25th  order  denominator,  as  presented  in  Appendix  A.  Many  of  these  roots 
occur  at  frequencies  higher  than  10  rad/sec  and  will  primarily  contribute  to  an  equivalent  time  de¬ 
lay.  Eight  of  the  denominator  roots  lie  in  the  frequency  range  of  .1  to  10  rad/sec,  making  it  dif¬ 
ficult  to  identify  the  dominant  modes  of  these  responses.  Approximate  pole-zero  cancellation  of 
roots  in  this  range  results  in  0/2nd  order  representations  for  both  transfer  functions.  Previous  ex¬ 
perience  leads  to  the  expectation  that  equivalent  system  models  will  be  easily  obtained  from  these 
representations.  A  good  sideslip  match  can  be  expected  while  the  bank  angle  response  should 
include  a  large  time  delay  due  to  the  additional  (uncancelled)  root  in  the  high  order  model. 

The  expected  equivalent  system  results  were,  in  fact,  obtained,  as  shown  in  table  VIII.  The 
approximate  forms  were  readily  matched  and  exhibit  relatively  low  mismatch  values.  The  roll  rate 
approximation  exhibits  a  large  time  delay  and  a  roll  mode  time  constant  different  from  that  ob¬ 
tained  from  simply  tracking  the  aircraft  root  in  the  high  order  system  representation.  The  equiva¬ 
lent  Dutch  roll  damping  is  similar  to  that  of  the  high  order  oscillatory  root  while  the  equivalent  fre¬ 
quency  is  somewhat  higher.  Simultaneous  roll  and  sideslip  angle  matching  (with  all  roots  in  the  fre¬ 
quency  range  of  .1  to  10  rad/sec  allowed  to  go  free)  resulted  in  excellent  matches  as  shown  in 
table  VIII.  Similar  trends  to  the  approximate  mode  were  obtained,  except  that  the  Dutch  roll  fre¬ 
quency  mode  closely  matched  that  of  the  high  order  system  oscillatory  root. 

Frequency  and  time  history  matches  for  the  simultaneous  match  results  are  presented  in  fig¬ 
ures  20  and  21.  Only  minor  differences  between  the  high  order  and  low  order  equivalent  system 
models  are  evident,  as  is  to  be  expected  from  the  low  (ss  4.5)  mismatch  values  obtained. 
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TABLE  VIII 


F-18  AIRPLANE  EQUIVALENT  TRANSFER  FUNCTIONS 
Roll  and  Sideslip  Angle  Response  to  Cockpit  Control  Inputs 
CR  IQ.OOOFT  .50  M 


'Fixed  parameter 
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Figure  21.  F - 1 8  Sideslip  Response  —  Simultaneous  Match 
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COMPARISON  WITH  SPECIFICATION  REQUIREMENTS 

The  lateral  directional  modal  parameters  resulting  from  the  equivalent  system  analysis  were 
compared  against  the  requirements  of  reference  (f)  to  determine  specification  compliance.  Consis¬ 
tent  trends  were  obtained  for  the  equivalent  system  parameters  resulting  from  the  approximate  and 
complete  formats  for  all  five  aircraft  analyzed  as  shown  in  figures  22,  23,  24,  and  25. 

There  are  differences  evident  in  the  modal  parameters  resulting  from  the  approximate  and 
complete  equivalent  system  formats.  In  both  the  cruise  and  power  approach  configurations,  the 
equivalent  Dutch  roll  damping  exhibited  the  largest  differences,  with  the  complete  format  results 
being  from  5  to  25%  higher  than  those  for  the  approximate  form.  The  complete  form  equivalent 
Dutch  roll  frequency  was  slightly  reduced  from  that  of  the  approximate  form,  while  roll  mode  time 
constant  showed  a  slight  increase  for  the  approximate  form.  The  differences  in  roll  mode  time  con¬ 
stant  were  only  significant  in  two  cases:  1)  the  A-6  airplane,  for  which  the  complete  form  roll 
mode  time  constant  was  consistently  0.2  seconds  or  more  greater  than  that  of  the  approximate 
form,  and  2)  the  A-7  airplane  at  0.3  Mach  number. 

All  conditions  analyzed  met  the  Category  A,  Level  1,  Dutch  roll  frequency  and  damping  re¬ 
quirements  of  MIL-F-8785C.  The  F-14  and  F-1 8  aircraft  met  the  configuration  CO  (Combat)  and 
G  A  (Ground  Attack)  requirement  that  damping  ratio  be  greater  than  0.4,  at  all  flight  conditions 
analyzed,  while  the  other  three  aircraft  only  meet  this  requirement  at  the  highest  Mach  numbers 
investigated. 

The  requirement  that  roll  mode  time  constant  be  less  than  1.0  second  was  met  by  all  five  air¬ 
planes,  at  all  conditions  analyzed,  except  one;  the  A-7  airplane  at  0.3  Mach  number.  This  configura¬ 
tion  exhibited  the  largest  discrepancy  between  the  approximate  and  complete  formats  and  was  also 
the  most  difficult  configuration  to  match  due  to  the  proximity  of  it's  higher  order  system  roots. 

Equivalent  time  delays  also  exhibited  consistent  trends.  The  complete  format  generally  re¬ 
sulted  in  lower  roll  and  sideslip  angle  time  delays  than  determined  for  the  approximate  forms.  The 
differences  were,  however,  small  and  of  little  consequence  in  the  resulting  time  history  responses  to 
control  inputs. 

The  identified  time  delays  met  the  Level  1  requirements  of  MIL-F-8785C  at  all  conditions 
analyzed  except  for  roll  angle  commands  in  the  A-7  and  F-18  aircraft.  In  the  case  of  the  A-7, 
presented  in  figure  26,  these  large  time  delays  arise  from  the  lateral  command  augmentation  roots 
which  are  not  explicitly  modelled  in  the  equivalent  system  format,  as  discussed  previously.  Flight 
experience  with  the  A-7  airplane  does  not  indicate  any  handling  quality  discrepancies  which  might 
be  expected  from  time  delays  of  this  magnitude. 

The  F-18  aircraft,  as  modelled  in  this  analysis,  utilizes  control  force  inputs  in  a  highly  complex 
digital  flight  control  system.  The  stick  force  filters  and  sensors  necessary  to  transmit  the  com¬ 
manded  force  to  the  flight  control  system  result  in  the  large  equivalent  time  delay  shown  in  fig¬ 
ure  22.  The  latest  revision  to  the  F-18  flight  control  system  does  not  use  control  force  inputs.  In¬ 
stead,  control  position  inputs  are  utilized,  eliminating  the  need  for  force  prefilters.  The  removal  of 
these  prefilters  has  reduced  the  time  delays  to  less  than  0.10  seconds,  which  fall  within  the  Level  1 
boundaries  as  defined  by  reference  (f). 

A  lateral  high  order  system  test  program  flown  on  the  NT-33  airplane  (reference  (q))  suggests 
that  roll  mode  time  constant  be  plotted  against  time  delay  to  determine  lateral  flying  qualities 
characteristics.  The  time  delay  parameter  in  reference  p  was  determined  graphically  and  is  slightly 
different  from  that  arising  from  the  equivalent  system  method  used  herein.  Therefore  the  configur¬ 
ations  of  reference  (q)  were  analyzed  by  the  present  methodology  and  the  boundaries  of  reference  (q) 
replotted  as  shown  in  figure  27  along  with  the  approximate  form  results  obtained  in  this  analysis. 

The  A-6,  F-14  and  S-3  aircraft  results  all  lie  within  the  Level  1  region,  while  the  A-7  and  F-18  re¬ 
sults  lie  in  the  Level  2—3  regions. 
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Figure  25.  Dutch  Roll  Mode  Results  -  Power  Approach  Configuration 


Figure  26.  A-7  Time  Delay  Results 
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Figure  27.  Time  Delay  vs.  Roll  Mode  Time  Constant  -  Approximate  Form  Results 
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CONCLUSIONS  AND  RECOMMENDATIONS 

Lateral  directional  equivalent  system  models  have  been  determined  for  five  tactical  Navy  air¬ 
craft  using  both  single  degree  of  freedom  roll  and  Dutch  roll  approximations  and  complete  three 
degree  of  freedom  formats  for  the  low  order  systems.  Mismatch  between  the  high  and  low  order 
systems  was  improved  by  using  the  complete  format.  However,  the  resulting  approximate  and  com¬ 
plete  modal  parameters  were  not,  with  the  possible  exception  of  Dutch  roll  damping  ratio,  signifi¬ 
cantly  different.  When  only  denominator  characteristics  are  desired,  the  single  degree  of  freedom 
roll  mode  and  Dutch  roll  approximate  formats  should  be  used  to  match  roll  rate  and  sideslip  angle, 
respectively.  These  approximate  models  will  provide  acceptable  equivalent  system  parameters  for 
comparison  against  the  MIL-SPEC  requirements. 

For  those  instances  in  which  numerator  characteristics  are  desired,  the  complete  three  degree 
of  freedom  equivalent  system  formats  must  be  used.  For  these  cases,  an  iterative  matching 
procedure,  in  which  numerator  and  denominator  roots  are  alternately  fixed  and  freed  has  been 
developed  to  optimize  the  matching  procedure. 

The  identified  equivalent  system  parameters  reflect  Level  1  flying  qualities  for  the  conditions 
analyzed  when  compared  against  the  requirements  of  MIL-F-8785C  with  the  exception  of  roll  angle 
time  delays  for  the  A-7  and  F-18  airplanes.  These  results  indicate  Level  2—3  equivalent  time  delays. 

It  is  recommended  that  these  frequency  response  matching  techniques  be  applied  to  other 
types  of  aircraft  to  determine  their  low-order  equivalent  systems.  Specific  aircraft  types  for  which 
these  techniques  would  be  applicable  include  VSTOL,  rotary  wing,  and  large  aircraft.  These  analyses 
should  include  both  the  longitudinal  and  lateral-directional  axes  and  would  provide  data  for 
examining  aircraft  flying  qualities  requirements. 
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APPENDIX  A 

AIRCRAFT  AND  CONTROL  SYSTEM  DESCRIPTIONS 


Data  describing  each  of  the  airplanes  and  their  respective  control  systems  were  obtained  from 
available  stability  and  control  reports,  references  (g)  thru  (I).  This  appendix  briefly  describes  the 
subject  airplanes  and  presents  a  block  diagram  of  their  respective  lateral  directional  control  system 
as  modelled  in  this  analysis. 
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S-3  -  The  S-3  airplane  is  a  twin  turbofan  powered,  land  and  carrier  based,  subsonic,  antisubmarine 
warfare  aircraft.  Lateral  directional  control  is  accomplished  via  a  mechanical  control  system  which 
operates  the  ailerons,  rudder,  and  spoilers.  Control  stick  dynamics  were  not  included  in  this  model 
of  the  S-3  aircraft.  The  aircraft's  basic  stability  is  augmented  through  the  feedback  of  both  roll 
rate  and  yaw  rate.  A  block  diagram  of  the  S-3's  lateral  directional  control  system,  as  modelled  in 
this  analysis,  is  presented  in  figure  A-1.  The  transfer  functions  representing  the  S-3  airplane  are 
presented  in  table  A-l. 
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A -6  -  The  A-6  airplane  is  a  twin  turbojet,  land  and  carrier  based,  subsonic,  all-weather  attack  air¬ 
craft.  A  block  diagram  of  the  A-6's  lateral  directional  control  system,  as  modelled  in  this  analysis,  is 
presented  in  figure  A-2.  The  control  stick  and  rudder  pedals  are  linked  directly  to  their  correspond¬ 
ing  surface  actuators  by  a  system  of  pushrods,  bellcranks,  and  cables.  Lateral  control  is  obtained 
through  the  use  of  flaperons  while  directional  control  is  provided  by  a  single  rudder.  The  basic 
stability  of  the  airplane  is  augmented  through  the  feedback  of  roll  rate  and  yaw  rate  to  obtain  the 
desired  response. 

The  transfer  functions  representing  the  A-6  airplane's  response  to  cockpit  control  deflections, 
obtained  via  the  Boeing  Computer  Services  program,  EASYS  (reference  r),  are  presented  in  table 
All. 
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A-6  AIRPLANE  TRANSFER  FUNCTIONS 
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A-7 
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A-7  -  The  A-7  airplane  is  a  single  place  turbo  fan  powered,  land  and  carrier  based,  light  attack  air¬ 
craft.  An  irreversible  mechanical  system  is  utilized  to  produce  lateral  directional  control  with  both 
stability  and  control  augmentation.  The  stability  augmentation  system  provides  roll  rate  and  yaw 
rate  feedback  signals  to  augment  the  aircraft's  basic  stability  characteristics.  The  command  augmen¬ 
tation  system  feeds  control  force  signals  forward  through  a  prefilter  as  a  means  of  increasing  the 
pilot's  commanded  input.  A  block  diagram  of  the  A-7's  lateral  directional  control  system,  as 
modelled  in  this  analysis,  is  presented  in  figure  A-3.  The  transfer  functions  representing  the  A-7 
airplane's  response  to  pilot  force  commands  are  presented  in  table  A-l  II. 
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F-14  -  The  F-14  airplane  is  a  twin  turbo-fan  powered,  land  and  carrier  based,  supersonic  fighter  air¬ 
craft.  Lateral  directional  control  is  accomplished  via  an  irreversible  mechanical  flight  control  sys¬ 
tem  which  transmits  cockpit  control  commands  to  a  differential  stabilizer,  spoilers,  and  rudders. 

The  airplane's  basic  stability  is  augmented  through  the  feedback  of  roll  rate,  yaw  rate,  and  lateral 
acceleration  to  obtain  the  desired  response.  A  block  diagram  of  the  F-14's  lateral  directional  con¬ 
trol  system,  as  modelled  in  this  analysis,  is  presented  in  figure  A-4.  The  transfer  functions  repre¬ 
senting  the  F-14's  response  to  cockpit  control  inputs  are  presented  in  table  A-IV. 
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Notes:  (1)  Direct  Lift  Control  (DLC)  ON 
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F-18  -  The  F-18  is  a  single  place,  turbo  powered,  land  and  carrier  based,  fighter  aircraft  controlled 
by  a  digital  flight  control  system.  Separate  flight  control  laws  are  provided  for  differeing  flight 
regimes.  Electrical  signals  are  generated  from  the  pilot's  control  force  inputs,  passed  to  the  com¬ 
puter,  modified  by  various  gains  and  shaping  networks,  and  finally  passed  to  the  rudder,  aileron, 
and  differential  tail.  Lateral  acceleration,  roll  rate,  and  yaw  rate  signals  are  also  input  to  the  com¬ 
puter,  where  they  are  shaped  and  gain  scheduled  before  being  summed  with  the  command  input 
signals.  The  cockpit  control  feel  system  dynamics  were  not  included  in  this  analysis.  A  simplified 
block  diagram  of  the  F-18's  lateral  directional  control  system,  as  modelled  in  this  analysis,  is  pre¬ 
sented  in  figure  A-5.  The  transfer  functions  representing  the  airplane's  response  to  pilot  commands 
are  presented  in  table  A-V. 
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Figure  A-5.  F-18  Lateral-Directional  Control  System 
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APPENDIX  B 

SUMMARY  OF  EQUIVALENT  SYSTEM  MODELS 


The  equivalent  system  results  for  the  S-3,  A -6,  A-7  and  F-14  configurations  investigated  are  sum¬ 
marized.  (The  single  F-18  configuration  investigated  is  presented  in  the  body  of  the  report.)  The 
identified  modal  parameters  and  quantitative  mismatch  values  are  presented  in  tables  B-1  through 
B-IV,  respectively.  Graphical  comparisons  of  high-  and  low-order  frequency  and  time  responses,  for 
those  configurations  not  previously  shown  in  the  Results  and  Discussion  section,  are  presented  in 
figures  B-1  through  B-1 8. 
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TABLE  B-l 

S-3  AIRCRAFT  EQUIVALENT  PARAMETER  RESULTS 
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70.7 

.554 

T01 

- 

- 

-  60.64* 

- 

- 

-1605.01* 

tp2 

- 

- 

.400 

- 

- 

.327 

TP3 

- 

_ 

.015* 

- 

- 

.007* 

lP 

- 

.013 

.034 

- 

.024 

.040 

Tr 

.312 

- 

.355 

.183 

- 

.248 

’’s 

- 

- 

166.69* 

- 

- 

119.6* 

fDR 

- 

.28 

.29 

- 

.47 

.52 

WDR 

- 

2.14 

2.08 

- 

3.71 

3.43 

18.2 

- 

1.8 

wm 

2.4 

M/j 

- 

14.8 

2.2 

- 

mm 

1.8 

‘Parameter  fixed  at  HOS  value 
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TABLE  B-ll 

A-6  AIRCRAFT  EQUIVALENT  PARAMETER  RESULTS 


CRUISE  CONFIGURATION  -  20,000  ft  ALTITUDE 

0.40  Mach 

0.72  Mach 

0.87  Mach 

APPROX 

0  P 

COMP 

0+0 

APPROX 

0  P 

COMP 

0  +  0 

APPROX 

0  P 

COMP 

0+0 

K0 

.670 

- 

.646 

2.559 

- 

2.04 

4.26 

- 

3.35 

- 

- 

.194 

- 

- 

.473 

- 

- 

.369 

“0 

— 

- 

1.56 

— 

- 

2.04 

— 

- 

2.80 

x<t> 

.0076 

- 

.0077 

.0144 

- 

0 

.022 

- 

.007 

kp 

- 

.0293 

.0003 

- 

.859 

1.004 

.0066 

% 

— 

- 

-131.58  * 

- 

- 

- 

- 

175.44  * 

Th 

— 

— 

.692 

— 

— 

- 

— 

.417 

Th 

- 

- 

.010  * 

- 

.007  * 

- 

- 

.006  * 

XP 

- 

.025 

.042 

- 

.036 

.055 

- 

.033 

.048 

a 

.461 

- 

.642 

.278 

- 

.873 

.166 

- 

.349 

Ts 

- 

- 

105.26  * 

- 

- 

- 

- 

144.09  * 

^DR 

- 

.251 

.299 

- 

.334 

.440 

- 

.366 

.420 

WDR 

- 

1.736 

1.71 

- 

3.063 

2.99 

- 

3.584 

3.46 

M. 

0 

145.4 

- 

121.8 

70.5 

- 

49.9 

41.6 

- 

31.9 

M0 

- 

m 

4.2 

- 

2.3 

8.6 

1.4 

2.8 

... 

*  F  ixed  Parameter 
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Figure  B-5.  A-6  Roll  Response  -  Simultaneous  Match,  .87  Mach 
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TABLE  B-lll 

A-7  AIRCRAFT  EQUIVALENT  PARAMETER  RESULTS 


CRUISE  CONFIGURATION  -  15,000  ft  ALTITUDE 

0.30  Mach 

0.60  Mach 

0.90  Mach 

APPROX 

0  0 

COMP 

0  +  0 

APPROX 

0  0 

COMP 

0  +  0 

APPROX 

0  0 

K0 

2.49 

- 

1.87 

7.9 

- 

5.97 

El 

- 

5.19 

?0 

- 

.51** 

- 

- 

.69** 

- 

.66** 

<^>0 

- 

1.64** 

- 

- 

2.56** 

- 

4.64** 

t0 

.267 

.247 

.237 

- 

.205 

.202 

- 

.166 

K0 

- 

.002 

.00003 

- 

.0061 

.00006 

- 

.010 

.00006 

T01 

- 

- 

-34.5* 

- 

- 

-1250.0* 

- 

- 

526.3* 

CM 

- 

- 

.72** 

- 

- 

.64** 

- 

- 

.51** 

t&3 

- 

- 

.015* 

- 

- 

.009* 

- 

- 

.006* 

*0 

- 

.025 

.032 

- 

.032 

.045 

- 

.038 

.052 

Tr 

.73 

- 

1.24 

.60 

- 

.59 

.50 

- 

.44 

B 

- 

- 

22.2* 

- 

- 

370.4* 

- 

- 

63.29* 

fDR 

- 

.29 

.25 

- 

.40 

.42 

- 

.45 

.50 

WQR 

- 

1.60 

1.83 

- 

2.29 

2.27 

- 

3.80 

3.69 

M0 

266.1 

- 

109.5 

84.8 

- 

31.2 

45.3 

- 

6.7 

M0 

75.8 

33.3 

- 

10.4 

■B 

- 

9.2 

2.2 

•Parameter  fixed  at  HOS  value 
••Parameter  iteratively  fixed  in  matching  process 


B-10 


EQUIVALENT  SYSTEM 


Figure  B-8.  A-7  Sideslip  Response  —  Simultaneous  Match,  .3  Mach 


F-14  AIRCRAFT  EQUIVALENT  PARAMETER  RESULTS 
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PA  CONFIGURATION 
-SEA  LEVEL 

121  KEAS  j 

1? 

8- 

.33 

.51 

.90 

.053 

o 

o 

-6.02 

1.58 

.008* 

00 

q 

.973 

-15.6* 

.36 

1.14 

80 

5.4 

X  <a 
O 

DC 

a. 

a. 

<0. 

1 

1 

1 

1 

*r*> 

1 

1 

1 

o 

1 

1 

n 

1.26 

1 

109.3 

.374 

1 

1 

.076 

l 

1 

1 

1 

1 

.575 

1 

_ 

1 

1 

8’86 

1 

CRUISE  CONFIGURATION  -  15.000  ft  ALTITUDE 

0.91  Mach 

COMP 

<P  +  0 

2.15 

.64 

4.37 

.047 

00 

8 

500.0* 

05 

O 

q 

.011* 

.064 

.549 

-3333.3* 

LO 

3.89 

4.1 

X  ca. 

O 

DC 

CL 

1 

1 

1 

1 

00 

00 

q 

1 

1 

1 

.049 

1 

1 

LV 

3.97 

i 

4.6 

2.74 

1 

1 

.072 

! 

1 

1 

1 

1 

.551 

1 

1 

12.6 

1 

0.795  Mach 

s'* 

8* 

2.14 

CM 

r- 

3.29 

00 

CO 

o 

.010 

* 

o 

o' 

o 

CM 

.671 

.013* 

.067 

.466 

-833.3* 

.64 

3.12 

80 

5.6 

APPROX 

P  P 

1 

i 

1 

1 

CM 

00 

q 

1 

1 

1 

.045 

1 

1 

.54 

3.49 

1 

8.2 

2.35 

i 

1 

.048 

i 

1 

1 

1 

1 

.476 

1 

1 

1 

3.6 

1 

0.715  Mach 

8 

2.00 

00 

r>. 

2.53 

.036 

o 

o 

-142.9* 

.846 

.015* 

.067 

.489 

-333.3* 

00 

CD 

2.48 

o 

5.8 

x  er. 

o  «*• 

DC 

CL 

a. 

<  CL 

1 

i 

1 

1 

.586 

1 

1 

1 

.041 

1 

1 

.544 

2.98 

1 

901 

o 

CM 

CM* 

i 

1 

CD 

o 

1 

1 

1 

1 

1 

ZLV 

1 

1 

1 

3.8 

1 

0.40  Mach 

CL  ~ 
5  oa- 

O  t 
a 

.64 

ZL' 

1.04 

.045 

900 

* 

00 

CO 

1.94 

.020* 

in 

o 

.701 

62.5* 

.59 

1.06 

ID 

APPROX 

P  0 

1 

1 

1 

1 

.267 

1 

1 

1 

.020 

1 

1 

.491 

1.515 

i 

38.0 

.683 

1 

1 

.054 

! 

1 

1 

1 

1 

.671 

1 

1 

1 

12.6 

i 

-©■ 

X 

3 

<53. 

* 

ca 

CM 

21 

CO 

21 

V- 

</) 

DC 

O 

•^5 

ac 

a 

3 

5 

ca 

5 

B-15 


Parameter  fixed  at  HOS  value 


.007)  <4.99)  (2.01  (18.41)  (66.06)  (20.0)  (20.0) 


581)  (7.17)  (20.0)  (27 JO)  1.69,4.321 
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